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Abstract 

We review the status of a quasi - degenerate doublets in nuclei, called pseu- 
dospin doublets, which were discovered about thirty years ago and the origins 
of which have remained a mystery, until recently. We show that pseudospin 
doublets originate from an SU(2) symmetry of the Dirac Hamiltonian which 
occurs when the scalar and vector potentials are opposite in sign but equal in 
magnitude. Furthermore, we survey the evidence that pseudospin symmetry 
is approximately conserved for a Dirac Hamiltonian with realistic scalar and 
vector potentials. We briefly discuss the relationship of pseudospin symmetry 
with chiral symmetry and the implications of pseudospin symmetry for the 
antinucleon spectrum in nuclei. 
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I. INTRODUCTION 



Dick Slansky and I met almost thirty years ago as assistant professors at Yale University. 
Although,for the most part, Dick's published research at that time was in hadron collisions 
IJ, he became interested in group theory and symmetry, particularly applied to elementary 
particles physics, during his stay at Yale. I had similar interest in group theory and sym- 
metry, but primarily applied to the many - body problem, particularly in nuclear physics. 
In the mid seventies Dick and I each accepted positions at Los Alamos National Laboratory 
and we continued to discuss symmetry in physics. At Los Alamos he published a number 
of papers in group theory including his very popular review on "Group Theory for Unified 
Model Building" [fj], and also published a book on affine Lie algebras ||. In 1989 he became 
Theoretical Division Leader at the Los Alamos National Laboratory, and, technically, my 
"boss". 

At the time we were at Yale I had heard about an observation that certain single - 
nucleon levels in spherical nuclei were found to clump into quasi - degenerate doublets ||,|5| . 
The term pseudospin doublets was applied to these states, although the quasi - degeneracies 
we assumed to be accidental. About two years ago I rexamined these pseudospin doublets, 
the origin of which remained a mystery, and I discovered that they are a consequence of a 
relativistic symmetry I discussed this revelation with Dick and he suggested that this 
symmetry may have some connection with chiral symmetry. In this paper I would like to 
discuss the progress that has been made in understanding pseudospin symmetry, including 
an enticing connection with chiral symmetry via QCD sum rules. 
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II. PSEUDOSPIN SYMMETRY 



The spherical shell model orbitals that were observed to be quasi - degenerate have non - 
relativistic quantum numbers (n r , £, j = £+1/2) and (n r — 1,£ + 2, j = £ + 3/2) where n r , £, 
and j are the single-nucleon radial, orbital, and total angular momentum quantum numbers, 
respectively . This doublet structure is expressed in terms of a "pseudo" orbital angular 
momentum £ = £ + 1, the average of the orbital angular momentum of the two states in 
doublet and "pseudo" spin, s = 1/2. For example, (n r si/2, (n r — l)d 3 / 2 ) will have £ = 1 , 
(n r p 3 / 2 , (n r — l)/5/2) will have £ = 2, etc. These doublets are almost degenerate with respect 
to pseudospin, since j = £ + s for the two states in the doublet; examples are shown in 
Figure 1. Pseudospin "symmetry" was shown to exist in deformed nuclei as well |||9| and 
has been used to explain features of deformed nuclei, including superdeformation |L0| and 



identical bands ||1T| , ^L2| . However, the origin of pseudospin symmetry remained a mystery 
and "no deeper understanding of the origin of these (approximate) degeneracies" existed 
I3| . A few years ago it was shown that relativistic mean field theories gave approximately 



the correct spin orbit splitting to produce the pseudospin doublets Hl4| . In this paper we 
shall review more recent developements that show that pseudospin symmetry is a relativistic 
symmetry P,|l5|,|7l| . 



III. SYMMETRIES OF THE DIRAC HAMILTONIAN 

The success of the shell model implies that nucleons move in a mean field produced by 
the interactions between the nucleons. Normally, it suffices to use the Schrodinger equation 
to describe the motion of the nucleons in this mean field. However, in order to understand 
the origin of pseudospin symmetry, we need to take into account the motion of the nucleons 
in a relativistic mean field and thus use the Dirac equation. The Dirac Hamiltonian, H, with 



3 



an external scalar, Vs, and vector, Vy, potentials is given by: 

H = a-p + P(m + V s ) + V v , (1) 



where we have set % — c — 1, a, (3 are the usual Dirac matrices | 16fl , m is the nucleon mass, 
and p is the three momentum. The Dirac Hamiltonian is invariant under an SU(2) algebra 
for two limits: Vg — Vy = constant and Vs + Vy = constant. For finite nuclei the constant 
will be zero since each potential will go to zero at large distances [|17]. (For infinite nuclear 
matter the constant could be non-zero.) 

The generators for the SU(2) algebra, Si, which commute with the Dirac Hamiltonian, 
[H , Si) — 0, for the case when Vs = Vy are given by (17 



_ a-psja-p (1-/3) A (l+_g) , , 



where = crj/2 are the usual spin generators and <ii the Pauli matrices. This reduces to 



where 



a* 



Si = U p §i Up = 23 P Pi - §i. (4) 



In (|) Up — — is the momentum-helicity unitary operator introduced in [14[] that accom 



plishes the transformation from the normal shell model space to the pseudo shell model 
space while preserving rotational, parity, time - reversal, and translational invariance. This 
symmetry limit leads to spin doublets, but we know that the spin orbit splitting in nuclei is 
very large so this limit is not applicable to nuclei. However we shall see that it is relevant 
for the antinucleon spectrum. 

We shall show that, in the limit of = — Vy, the conserved symmetry is pseudospin 
symmetry. The generators for the SU(2) algebra, Si, which commute with the Dirac Hamil- 
tonian, [ H , Si ] — 0, for the case when Vs = —Vy are given by |L7 



~ = a-psjd-p (1 + P) (1 - P) 

i p2 2 Si 2 U 

This reduces to 

For V$ = —Vy, the eigenf unctions of the Dirac Hamiltonian, H^/ T jx = £ t ^t,[i are doublets 
(S — 1/2, fx — ±1/2) with respect to the SU(2) generators S{ of Eq. (0) 



S± %,fi = ^(1/2 =F /2)(3/2 ± /2) * T)A± i , (Vs = -V v ) (7) 

where 5*± = S x ± z^y. The eigenvalue r refer to the other necessary quantum numbers. 

The fact that the the pseudo - spin generators (|) have only the spin operator Sj operating 
on the lower component of the Dirac wave function has the consequence that the spatial 
wavefunctions for the two states in the pseudospin doublet are identical in the limit of 
Vs — —Vy to within an overall phase. 

This symmetry for Vs — — Vy is general and applies to deformed nuclei as well as spherical 
nuclei. In the case for which the potentials are spherically symmetric, the Dirac Hamiltonian 
has an additional invariant SU(2) algebra; namely, the pseudo - orbital angular momentum, 









(8) 



where l{ = U p £i U p , li — r x p. In this limit, the Dirac wave functions are eigenfunctions 
of the Casimir operator of this algebra, L ■ L l^^,^.) = + l)|^ T ijm-)> wri ere we have 
used a coupled basis, and j is the eigenvalue of the total angular momentum operator 

Ji = Li + Si, J ■ J\^ Tt ij,m ) = j(J + l)l*T/,i,m )> anci m i is tne eigenvalue of J z . Thus 
pseudo-orbital angular momentum as well as pseudospin are conserved in the spherical limit 
and Vs = —Vy. From (Rt), we see that the lower component wave function will have spherical 



harmonic of rank t coupled to spin to give total angular momentum j. Since a ■ p conserves 
the total angular momentum but p changes the orbital angular momentum by one unit 
because of parity conservation, the upper component also has total angular momentum j, 
but orbital angular momentum £ = I ± 1. If j = £+1/2, then it follows that I = I + 1, 
whereas if j = i — 1/2, then i = £ — 1. This agrees with the pseudospin doublets originally 
observed |4],§] and discussed at the beginning of this paper. 

For axially symmetric deformed nuclei, there is a U(l) generator corresponding to the 
pseudo-orbital angular momentum projection along the symmetry axis which is conserved 
in addition to the pseudospin for Vs = —Vy, 



value for the upper and lower components since a • p conserves the total angular momentum 
projection. Thus Q = A ± 1/2, corresponding exactly to the quantum numbers of the 
pseudospin doublets for axially deformed nuclei discussed in 0|§. 

However, the exact symmetry limit can not be realized in nuclei, because, if Vs = —Vy, 
there are no Dirac bound valence states and hence nuclei can not exist. However, we now 
show that pseudospin symmetry is approximately realized for Vs ~ — Vy. 

IV. FREE NUCLEONS 

For a free Dirac nucleon, Vs = Vy = 0, and hence both symmetries are valid. The free 
nucleon wavefunction is 




(9) 



where A = U p A U p . In this case the Dirac wave functions are eigenfunctions of A, A \^f r ^ n ) = 
M^tAu)^ where Q is total angular momentum projection, Q = A + fl, which has the same 




(10) 
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where are the Dirac spinors Ui/ 2 = (q),W-i/2 = (°), E p = ^p 2 + m 2 , and M is the 
normalization. This wavefunction has definite helicity and hence is an eigenfunction of the 
helicity operator, 

U p V Pttt = 2nV Ptlt . (11) 

The wavefunctions with \i = ±1/2 are doublets for Si and, since S% — U p SiU p ^ the wave- 
functions ^ p ,p, = Upfyp jlx = 2fi^/ PtP are doublets for Si, the pseudospin. 



V. REALISTIC MEAN FIELDS 



A near equality in the magnitude of mean fields, V$ ~ — Vy, is a universal feature of the 
relativistic mean field approximation (RMA) of relativistic field theories with interacting 
nucleons and mesons |18] and relativistic theories with nucleons interacting via zero range 



interactions ||19|| , as well as a consequence of QCD sum rules ||20|| . We shall discuss QCD 
sum rules in the next section. 

Recently realistic relativistic mean fields were shown to exhibit approximate pseudospin 



symmetry in both the energy spectra and wave functions ||T5| , pT| , |22| . In Fig. 2 we show the 
energy splittings between pseudospin doublets normalised by H + 1 as a function of the 
average binding energy (e) = (e| +1 / 2 + e £-i/2)/2- We see that the energy splitting for the 
same pseudo - orbital angular momentum decreases as the radial quantum increases; that 
is, as the binding energy decreases. Also for the same binding energy, the energy splitting 
increases as the pseudo - orbital angular momentum increases. These features follow from 
the square well potential 0. In Table 1 we tabulate some of the calculated energy splittings 
compared to the measured splittings. What is interesting is that the measured splittings are 
smaller than those calculated by relativistic mean field theory indicating that pseudospin 
symmetry breaking is overestimated by the relativistic mean field approximation. 
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Pseudospin doublets will manifest themselves for deformed potentials as well when Vs ~ 
—Vy. In Figure 3 the single particle (s.p.) energies of the doublets are plotted versus the 
deformation. We see that for each value of pseudo asymptotic quantum [A r n 3 A] numbers 
there is a quasi-degenerate pseudospin doublet with ft — A ± 1/2. 

As mentioned in the last section the relativistic SU(2) pseudo-spin symmetry im- 
plies that the spatial wavefunction for the lower component of the Dirac wave- 
functions will be equal in shape and magnitude for the two states in the dou- 
blet. For spherical nuclei the Dirac wavefunction for the two states in the dou- 
blet are ^ T j=i + \/2,m = (5ri+l/2K+lX]m <+1 ^)«/ T ,f+l/2K'X]m f+1 ^)) ^T,j=£-l/2,m = 

(^r,£-i/2[ Y £-iX]^~ 1/2) ^/ T ,£-i/ 2 [ y £X]^" 1/2) ) where g,f are the radial wave functions, Y e 
are the spherical harmonics, % is a two-component Pauli spinor, and [. . .] means coupled 
to angular momentum j. For a square well potential, the overall phase between the two 
amplitudes will be a minus sign |J so we expect that, in the symmetry limit for realistic po- 
tentials, f T i +1 / 2 ( r ) = ~f T i-i/2( r )- m Figure 4 we see that, for realistic zero range potentials, 



fr,e+l/2( r ) ~ -/r,£-l/ 2 ( r ) HI" 

These results are also valid for the relativistic mean field approximation to a nuclear field 



theory with meson exchanges [21]. In Figure 5a) the pseudospin doublets in the vicinity of 
the Fermi surface for neutrons and protons are shown. The upper (g) and lower components 
(f ) of the pseudospin doublets are also shown in Figure 5b,c,d. While the upper components 
are very different with different nodal stucture, the lower components are almost identical. 
We also note that the lower components are small with respect to the upper components, 
which is consistent with the non-relativistic shell model. 
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VI. QCD SUM RULES 



Applying QCD sum rules in nuclear matter, the scalar and vector self-energies were 



determined to be [20 



4n 2 a N p N 
m z m q 

K - ^SL, (12) 

where is the nuclear density, and m q the quark mass, a at is the sigma term which arises 
from the breaking of chiral symmetry PBfl . The ratio then becomes 

5i = — — . (13) 
8m q 

For reasonable values of <jn and quark masses, this ratio is close to -1. The implication of 
these results is that chiral symmetry breaking is responsible for the scalar field being approx- 
imately equal in magnitude to the vector field, thereby producing pseudospin symmetry. 

VII. ANTINUCLEON SPECTRUM 

The antinucleon states are obtained by charge conjugation, C, applied to the negative 
energy eigenstates of the Dirac Hamiltonian [[TjJ . This leads to a spectrum which has quasi 



- degenerate spin doublets, not pseudospin doublets. This follows from the fact that, under 
charge conjugation, 

c,| " c= (o1H" <i4) 

Thus in fllil) the spin operator Sj operates on the upper component and hence the spatial 
wavefunctions for the upper components of the states in spin doublet will be very similar. 
Likewise for spherical nuclei, the pseudo-orbital angular momentum goes into the orbital 
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angular momentum, and for axially deformed nuclei, pseudo-orbital projection goes into 
orbital projection along the body fixed z axis. 

This symmetry in the antinucleon spectrum also follows from the fact that the antinu- 
cleon potentials are Vs = CWs C = Vs, and Vy = CWy C = —Vy. Thus Vs ~ Vy and the 
symmetry of the Dirac Hamiltonian generated by (|3]) applies and spin doublets are produced 



in the antinucleon spectrum [1/ 



VIII. SUMMARY 

We have shown that pseudospin symmetry is a broken SU(2) symmetry of the Dirac 
Hamilonian which describes the motion of nucleons in realistic scalar and vector mean field 
potentials, Vg ~ — Vy. This symmetry predicts that the spatial wavefunctions of the lower 
components for states in the doublet will be very similar in shape and size and this has been 
substantiated by relativistic mean field approximations of relativistic nuclear field theories 
and relativistic nuclear Lagrangians with zero range interactions. This symmetry has been 
linked via QCD sum rules to chiral symmetry breaking in nuclei. Finally, the antinucleon 
spectrum is shown to have a spin symmetry rather than a pseudospin symmetry. 

Future applications of pseudospin symmetry will involve the testing of the wavefunction 
through the relationships between transitions within pseudospin doublets which follow from 
pseudospin symmetry ||24|| . 
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TABLES 

TABLE I. 208 Pb pseudospin doublet energy splittings 

1 ps doublets e l+1/2 - e< -_ 1/2 (RMA) |TJ e i+1/2 - q_ 1/2 (EXP) 
neutrons 

4 0h 9/2 - lf 7/2 2.575 1.073 

2 l/ 5/2 - 2p 3/2 0.697 -0.328 

protons 

3 0g 9/2 - ld 5/2 4.333 1.791 
1 ld 3/2 -2s 1/2 1.247 0.351 
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FIGURES 

FIG. 1. Examples of pseudo-spin doublets in the 208 Pb region. n r is the radial quantum number 
of the state, t is the orbital angular momentum, j the total angular momentum. 

FIG. 2. Pseudospin doublet energy splittings normalised by 2£ + 1 as a function of the average 
binding energy (e). n r is the radial quantum number of the state with the lower orbital quantum 
number 

FIG. 3. Single - particle (s.p.) energies for neutron pseudospin partners as a function of defor- 
mation 

FIG. 4. 208 Pb lower component wavefunctions f T £_i/ 2 ( T ) (dash line), -/ r £ +1 / 2 ( r ) (dot - dash 
line) for the (2s 1 / 2 , W3/2) pseudospin doublet as a function of the radius r 

FIG. 5. 208 Pb: a) energy spectrum; b,c,d) upper (g) and lower (f) components of the Dirac 
wavefunction for pseudospin doublets as a function of the radius r 



14 



Examples of Pseudospin Doublets 

(n r if, (n r -l)£ + 2j +1) 

j = I ± s , s = 1/2 ; 
£ pseudo-orbital angular momentum, s pseudo-spin 



3s 



1/2 



2d 



! j 

( 1 )3/2 , 1/2 



3/2 



3p, 
2f 



3/2 



( 2 ) 5/2 , 3/2 



5/2 



2d 



5/2 



ig 



(3) 



7/2 , 5/2 



7/2 



2f 



7/2 



1h 



(4)9/2, 7/2 



9/2 
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